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Abstract 

We present a straightforward method for obtaining exact classical and quantum molecular 
Hamiltonians in terms of arbitrary coordinates. As compared to other approaches the resulting 
expression are rather compact, the physical meaning of each quantity is quite transparent, 
and in some cases the calculation effort will be greatly reduced. We also investigate systems 
with constrains to find the suggested method to be applicable in contrast to most conventional 
approaches to kinetic energy operators which cannot directly be applied to constrained systems. 
Two examples are discussed in detail, the monohydrated hydroxide anion and the protonated 
ammonia dimer. 
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I. INTRODUCTION 



Theoretical modeling of moelcular reaction dynamics requires to have at hand a Hamil- 
tonian. For few-atom systems this may involve defining the potential energy surfaces 
(PES) explicitly in terms of properly chosen coordinates such as bond length and an- 
gles or distance vectors. In the latter case, for instance, it has been shown that PES 
in full-dimensionality including proper account of exchange symmetry can be generated 
for systems like the five-atomic HsO^ [1] or the seven-atomic HsO^ [2]. With increas- 
ing dimensionality this approach is no longer feasible due to the exponential scaling of 
the computational effort for explicit calculation of the PES on a numerical grid. Here, 
reaction surface concepts have been proven to be a versatile compromise, mapping out 
that PES region which is relevant for the considered reaction by combining suitable large 
a,npUtude .e.t.o coordinates witK o.tKogo.a. Mnnon.c .0*10.^0,0,1^3. 

Besides the PES the kinetic energy operator (KEO) needs to be known in terms of the 
chosen coordinate system. In principle, kinetic energy quantisation had been discussed by 
Podolsky in terms of intricate tensor analysis right after the foundation of the Schrodinger 
equation had been laid Later studies are to some extent based on the same mathemat- 
ical techniques [o], [l^. Approximate KEOs using normal mode coordinates are 



widely adopted which go back to contributions by Eckart Wilson fl2'|, and 
Watson 1 131]. In particular, the Eckart equations ll|] enable one to determine 
an orientation of molecule fixed axes suitable for using normal modes. The 
details of the derivation of a rovibrational Hamiltonian can be found in Ref. 



12[. A significant simplification of this rovibrational Hamiltonian in terms of 



normal coordinates has been achieved by Watson [13, 14]. In fact the Watson 
Hamiltonian served as a starting point for many later investigations and is 
implemented, e.g., in the MULTIMODE program of Bowman, Carter, and 
coworkers [15]. In principle normal mode coordinates provide a very compact 
representation of the PES which is taylored, however, to the stationary point 
for which they have been defined. For large amplitude motions away from 
this stationary point the above mentioned reaction surface approach or the 
use of general curvilinear coordinates are more suitable. A detailed discus- 
sion of advantages and disadvantages of various coordinates for rovibrational 
Hamiltonians has been presented, e.g., in Ref. |16l ]. 

In recent years Gatti and coworkers have developed a general scheme for obtaining a 



KEO using vector parameterization, i.e., the non-trivial coordinates are expressed with 



— 1 real vectors for system of atoms 



17 



18 



convenient for a full dimensional description 20||. However, often one needs a reduced 



191]. This scheme is quite successful and 



21 



Generally 



dimensional description leading to a problem with imposed constraints 
speaking in this case the coordinates can not be expressed as components of real vectors. 
Finally, we should point out that a numerical alternative to the analytical determination 



of the KEO has been suggested in Ref. 



23|. 



In the present work we introduce a simple physical method of kinetic energy quanti- 
sation which directly quantises the classical Lagrange/Hamilton functions. The formal 
theory is presented in Section [ITl Here it is also shown that our method can be directly 
applied to constrained systems. In Section IIIII the approach is applied to two molecular 
systems, that is, H3O2 which is treated in full dimensionality and N2H^ for which certain 
constraints are applied. 



II. GENERAL THEORY FOR OBTAINING A KEO 



A. Formal Development 

In the following we focus on the most common cases where the kinetic energy includes 
purely quadratic terms and the potential energy does not depend on the velocity, although 
a generalisation would be straightforward. Further we adopt the laboratory reference 
frame (LRF). The classical Lagrange/Hamiltonian in terms of arbitrary variables can be 
written as 

L(Q, Q, t) = T{Q) - V{Q) = ^Q^MQ - V{Q) 

HiQ, P) = T{P) + V{Q) = iptM-ip + y(Q), (1) 

where Q and P are single column vectors of the generalised coordinates and correspond- 
ing conjugate momenta, respectively. The generalised momentum vector is defined as 
P = dL/dQ = ISAQ. The generalised mass matrix, M, in general is a function of the co- 
ordinates. We assume that it is defined such as to be Hermitian and "f" means Hermitian 
conjugate which is equivalent to the transpose in classical mechanics. In order to obtain 
the quantum mechanical Hamiltonian operator, the major task is to derive the quantum 
KEO in coordinate representation. 

Since the choice of coordinates is quite arbitrary in Eq. ([1]), we start with Cartesian ones 
{X,Px}', in the following all the subscripts "x" are associated with Cartesian coordinates. 



The operator for each component of the momentum (-Px)i — {■^x)j = —ihd/dXj is 
Hermitian. For non-Cartesian coordinates we define the momentum operator in the same 
way: 



_d_ 



(2) 



To obtain the quantum KEO in terms of non-Cartesian coordinates, we must exploit the 
coordinate transform relations, Q = Q{X). The following relations can be derived easily 



Q = 




dX^ 






P = 




\dQ) 


dQ 


rax 


dX 


~ [dQ 



Px 

(3) 

where the elements of the transformation matrices are defined as {dQ / dX)^ - = dQi/dXj 
and {dX /dQ)^j = dXi/dQj. The quantum KEO in terms of non-Cartesian coordinates 
can be obtained by the simple coordinate transform from Cartesian coordinates 



T 



X 



dX' 



dQ 
dX 



(4) 



where Mx is the diagonal Cartesian mass matrix consisting of real mass associated with 
each Cartesian coordinate. 

The momentum operator vector P can be replaced by the derivative operator P = 
—ihd/dQ. However, the expression for the Hermitian conjugate of momentum operator 
(HCMO) vector P^ needs to be established yet. Since one has the relation between P 
and the Cartesian one Px it is not difficult to directly apply the definition of Hermitian 
conjugation to Eq. ([3]) 




(5) 



where the notation [A, ]3] means the commutator of two operators. Eq. shows the 
non-Hermiticity of generahsed momenta associated with non-Cartesian coordinates. Addi- 
tional terms appear due to the non-commutabihty of the momenta and the transformation 
matrix. Using the basic commutator [Xj, {Px)k] = i^Sjk, Eq. ([5]) can be simphfied to (see 
also Ref. [2^) 



Here, the superscript ° means a differential operator inside the bracket can not operate 
on functions outside, in other words, the result is just normal function of coordinates. 

Eq. ([6]) clearly shows the relation between a momentum operator and its Hermitian con- 
jugate. Apparently, any set of functions of coordinates, {fj{Qj)}, can be multiphed from 
the left to the derivative operator in Eq. ([2]) yielding different schemes of momenta quanti- 
sation, i.e., (-P)j = —ih- fjiQj) -d/dQj constitutes also an acceptable scheme. Since these 
momentum operators are in general non-Hermitian, we cannot simply set fj{Qj) = 1. For 
different schemes of momenta quantisation, the HCMOs will vary correspondingly to keep 
the KEO invariant. But, this enables us to optimise momenta quantisation schemes. In 
other words with proper functions one can obtain desired forms of momentum operators, 
e.g., symmetric forms. 

With Eq. (jl]) and Eq. we get the general scheme of constructing KEOs in terms of 
arbitrary coordinates. However, the final structure seems to be complicated and actually 
it can be simplified. To this end we compare Eq. (jl]) and Eq. ([1]) and take into account 
the invariance of the classical kinetic energy 

T = ^Q^MQ 

Since Q can be any vector we must have the relation M = [dX / dQ)'^ ^i^dX / dQ. 
Calculating the inverse of both sides we immediately see that the complicated central 
part (9Q/aX)Mx^ (dQ/dX)^ in Eq. (g]) is exactly M^^ in Eq. (P. The final KEO is 
therefore simplified as 

f = -P^M-^P. (8) 

This means we can directly exploit the result of Eq. ([1]) and quantise the generalised 
momenta without the knowledge of the Cartesian kinetic energy. All we need are the 
generalised mass matrix M and the definition of coordinates which both appear in Eq. ([1]). 



Therefore we can identify a concise and physically transparent scheme to construct the 
KEO: 



1. Get the classical kinetic energy and make sure the mass matrix is symmetric, so 
that one has the same structure as Eq. ([1]), T{Q) = ^Q^lSAQ. 

2. Exploit Eq. or Eq. (jH]) to express HCMOs in terms of P and some functions 
of Q. 

3. Calculate the inverse matrix of M hence the formal quantum KEO reads T = 
Ip^M^^P. 

4. Replace P by —ihd/dQ for obtaining the coordinate representation. 

The first step is quite familiar and one can choose arbitrary coordinates to get the classical 
kinetic energy. Actually we will introduce a useful partition method in Section Hi CI which 
will significantly simplify this issue. The remaining three steps are quite straightforward 
to follow. Apparently, the major effort is in the second step, i.e. determining the HCMOs 
according to Eq. ([6]). However, for the familiar spherical coordinates the HCMOs are 
well-known which leads to a great simplification as shown in Section HID! 

Apart from the KEO the volume element of integration is also of importance. Since we 
start from Cartesian coordinates and all that we have done is a coordinate transformation, 
the Euclidean normalisation remains correct: 

dr = drx = l[ dX, = |Det \ J] dQ,. (9) 

B. Systems with Constraints 

Reduced dimensional descriptions are always necessary for large systems. In the follow- 
ing we will give the general description for systems with constraints. Considering a system 
with some active coordinates Qi and some frozen coordinates Qq, the full dimensional 
coordinates and corresponding conjugate momenta read 





The constraining conditions are given by Qq = 0. Thus we can obtain the constrained 
classical kinetic energy 

T = ^Q^MQ 

1 / . t . t \ Mil Mio 

2 V " ^ \^ Moi Moo 

= ^qIMuQi, (11) 

where Mjj are corresponding sub-matrices. 

To get the quantum KEO we have to rewrite the constraint conditions in terms of 
momenta. According to the definition of momenta it is not difficult to find the following 
relation 

-1 





(12) 



, A B , 

where I I is the inverse of matrix M. Thus one can rewrite the constraint condi- 

C D 

tions as CPi + T)Po = or Pq = — D~^CPi. After quantization we get the constraint 
relation for the corresponding quantum operators, i.e., Pq = — D^^CPi. Based on this 
observation we can obtain the quantum KEO 

f = -P^M-^P 
2 

if pt ptWMii Mio 
-^"^^"^^l^ Moi Moo 




2 

1{pI (-D-iCPi)t 
ipl (A - BD-^C) Pi 

^plMfi^Pi. (13) 



The final result is quite compact. Comparing Eq. ffTTj) with Eq. (1131) we observe that 
the same procedure mentioned in the full dimensional case can be followed provided we 
only consider the active coordinates and completely ignore the frozen ones from the very 



beginning when we generate the classical Lagrangian. This is an attractive point since 
when we need to deal with constrained systems the present scheme only requires the 
classical Lagrangian which can be obtained by traditional methods. 

C. A General Method for Partitioning the Classical Kinetic Energy 

So far we have given the general theory of kinetic energy quantization starting from 
the classical Lagrangian. In the following we will introduce a method for obtaining the 
classical kinetic energy with a partitioning technique which will greatly simplify the prob- 
lem in most cases. In general it is quite convenient to divide a large system into small 
subsystems especially when a subsystems has certain symmetry. If we arbitrarily divide 
a system into parts the kinetic energy is a sum of the subsystems. According to the 
Konig theorem, 

N N 

T = Y.Tf = Ttc) + Y.Tf\ (14) 

i=l i=l 

where is the kinetic energy of the part A with respect to the reference frame defined 
by B. Here is the laboratory reference frame and (C) is the centre of mass reference 
frame. 

In the special case when N = 2 the following relation can be obtained 

N 

J2 tP = tP + tP = Tp'^ + Tp^ + (15) 

i=l 

where (Cj) is the centre of mass of the ith part and T^-^^^ is the kinetic energy of the 
centre of mass of the first part with respect to the centre of mass of the second part. If 
the two parts are both mass points Eq. f|T5l) can be simplified as follows 

rp{C) . rp{C) _ rp{C2) _ rp(Cl) I ^ a\ 

which is a familiar result from two-body mechanics. By exploiting Eq. f|T5l) repeatedly 
we can easily express the total kinetic energy in terms of kinetic energies of subsystems 
which are much easier to obtain as will be seen in the following. 

D. Hermitian Conjugates of Momentum Operators 



As mentioned in Section fllAl the most tedious task for generating a KEO is to derive 
the expressions for the HCMOs. As an example let us consider a system described by the 



3D spherical coordinates {R, 6, (f}. The coordinate transformation between 3D spherical 
coordinates and 3D Cartesian ones {x, y, z} is defined as 



X = R sin 6 cos (f 
y = i? sin 6' sin (y9 

z = RcosO. (17) 

Here, we only give the corresponding HCMOs leading to the well-known KEO (for details, 
see Appendix IV AI) 

1 d 

PI = Pe — ih cot 9 = —ih -— sin 6 

sm 9 o9 

d 

^ ^ dip 

T = ^— P^- t: r^sin^— - 



2mR^dR dR 2mR'^ sin 9 09 09 2mR'^ sin^ 9 O'p'^' ^ ' 

So far our considerations have been concerned with the LRF. However, in order to 
describe a molecule in terms of its natural motions, e.g., bond lengths and bond angles, 
it is usually more convenient to use one or more molecular reference frames (MRFs). In 
the following we will discuss this point and derive results for HCMOs associated with 
MRF spherical coordinates. Consider the LRF and a MRF defined by sets of unit vectors 
{ca;, e^^, e^} and {e^;/, Cj,/, Cz'}, respectively. The relation between the LRF and the MRF 
is just an orthogonal transformation characterised by the three Euler angles {■(9, 0, x} 

= U,(0)U^(t9)U,(x)e«, (19) 

where a = x,y,z and Uq is a rotation around Cq. (The expression for the rotational 
transformation matrix and more details on the following derivations can be found in 
the Appendix |V3) 

Now let us consider a vector Rj characterised by three spherical coordinates {Rj, 9j, 
(Pj} in the MRF. According to Appendix |VAl Eq. (ge]), we can express Rj as 

R, = i?,U,(0)U,(^?)U,(x)U,(^,)U,(0,)e,. (20) 

This is a vector equation and we can obtain the Cartesian coordinates of Rj in LRF by 
projecting the equation onto each LRF axis. Based on Eq. (!20|) we can exploit Eq. ([6]) to 
derive the expressions for HCMOs associated with MRF coordinates. The final results 



are quite concise and they actually have the same form as those associated with LRF 
spherical coordinates. That is to say, that for the momentum operators associated with 
Rj {Rj, 9j, ipj) the following relations are valid 





= Pn. 


2ih 
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— ih cot 9j 
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sin 6j 89 j ^ 






= P^j 


dip, 




(21) 



In Appendix IV Al we also confirm that the momentum operators associated with MRF 
Cartesian coordinates are Hermitian (cf. Eq. (1631) ). This result is quite important as it 
shows that one does not need to take the effort to derive the expressions for HCMOs 
provided that one uses spherical coordinates, Cartesian coordinates or combinations of 
both, no matter whether they are defined in the LRF or MRFs. 



E. Angular Momentum and Rotation 

For a system with atoms we have 3A^ degrees of freedom (DOFs). They are normally 
divided into three translational DOFs, three rotational and 3A^ — 6 vibrational ones. The 
three translational DOFs can be separated while the remaining 3A^ — 3 DOFs are coupled. 
In general rotational excitation energies are quite small as compared with vibrational ones 
suggesting an approximate separation. For the following discussion it will be convenient 
to express the KEO as sum of rotational (three angles) and vibrational (3A^ — 6 other 
coordinates) parts as well as their coupling. Specifically, we exploit the 3D vector Q^^^ 
defined as 

Qlt = ( ^ X ) (22) 
containing the three Euler angles which connect the LRF and MRF according to Eq. (1191) . 
The kinetic energy and total angular momentum J are defined as 

N 

2T = ^miklRi 

i=l 
N 

J = m,Ri X Ri . (23) 

i=l 

With Qrot being the angular velocity of MRF the velocities can be re-expressed as 

R, = xR, + i?:, (24) 



where Ft!^ is the velocity of Ri measured in the MRF. With the help of Eq. fl2^ and some 
vector algebra we can rewrite the kinetic energy and angular momentum as 

N 



1=1 

N 

J = Y.m, [r^Q^,, - {qIMR. + {R. X . (25) 



Using Eqs. fl25|) it is straightforward to derive the following relation 

dT 

which is exactly the definition of the generalised momentum vector Prot associated with 
the three Euler angles. 

Now we can draw the following important conclusion. If the set of coordinates contains 
the three Euler angles {"d, 0, x}, the total angular momentum vector is just the generalised 
momentum vector associated with the three Euler angles 







dT 
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= = 






J<f, 
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dT 


Jx 


= Px = 


dx 
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= e^Pg 


+ e 



, + (27) 

However, since the direction of angular velocities {e^, e^, e^} are complicated, it is better 
to transform the expressions to the Cartesian components in the LRF [l^, i.e. 

= sin xPi) — CSC -i? cos xP(i> + cot cos xPx 
Jy = cos xP-& + CSC {} sin xP(f) — cot {} sin xPx 

Jz = Px- (28) 

Using Eq. (128!) it is straightforward to express the KEO in terms of total angular mo- 
mentum and generalised momenta associated with vibrational DOFs. In other words, 
starting from Eq. (1251) the KEO can be readily expressed as contributions of rotational 
part, vibrational part and their coupling. 

If the focus is on the vibrational spectrum, there are several arguments for neglecting 
rotational excitations. There is usually a clear time scale separation, that is, rotational 
excitation energies are much smaller than vibrational ones. If we are not aiming at high 



resolution spectroscopy, rotational excitation will show up merely as a broadening of the 
vibrational bands and the band shifts of vibrational transitions due to the coupling can 
be considered as being much smaller than the accuracy achievable by the quantum chem- 
istry and quantum dynamics methods for obtaining the spectra. Therefore the following 
considerations will assume that rotational motion can be neglected, that is, we will set 
J = OT P^ = P^ = P^ = 0. 



III. APPLICATIONS 



In the following we will give two applications of the present approach. First, we 
will consider the full-dimensional KEO for H3O2 employing polyspherical coordinates 
as introduced in more detail in Appendix IV B[ Second, we consider N2H7 under the 
assumption of certain constraints and using bond length and angle coordinates. 



A. Full Dimensional KEO for HsO^ 

In this section we will give the full dimensional KEO for the monohydrated hydroxide 
ion, H3O2 , which had been used in the study of different isotopomers in Ref. 25|]. The 
coordinates are defined in the MRF given by Eq. ( |T9i) . The four Jacobi vectors shown in 
Fig.[T] are adopted. The three Euler angles are chosen in the same way as introduced in 
Section rVBl i.e., the MRF spherical coordinates -R4 and Ri are (i?4,0,0) and (i?i,6'i,0), 
respectively. The other two vectors can be characterised by their spherical coordinates 
in the MRF {Rj, 9j, ipj) (j = 2, 3). We can calculate the G matrix according to Eq. (jTHD 
and the required matrix elements are shown in Eq. (1781) . Therefore, we can directly write 
down the exact 9D KEO for total angular momentum J = as 2Tvib = -PlibGvib-Pvib 
with the HCMOs given in Section IIIDi However, this KEO will contain a large number 
of terms describing the angular momentum coupling between the shared Hydrogen and 
the (0H)2 fragments. In order to simplify the KEO we will derive an approximate KEO 
which has the advantage that these couplings do not appear. 

The basic idea is to introduce a new MRF and express -R3 with Cartesian coordi- 
nates {x,y,z). The remaining six coordinates are spherical coordinates in the old MRF 
(-Ri, i?2, Gi, O21 <^ = v^2)- The new MRF is associated with the old MRF by a rotation 
of an angle rjip around the (the 2;-axis). The spherical coordinates for in the new 
MRF are (i?3, 6*3, ips) with ip^ = ips — rjip, where i] = /ii/ (yUi + /i2) is defined to minimise 



the Coriolis type couplings involving the central Hydrogen. Here, /ii (/i2) is the reduced 
mass associated with the Jacobi vector Ri (-R2)- After transforming the spherical coor- 
dinates (-R3, 6^3, ips) to the Cartesian ones (x, y, z) we can obtain the KEO for the central 
Hydrogen 

T =-^(^ + ^ + ^\ (29) 

where /is is the reduced mass associated with the Jacobi vector R3 and Coriolis type 
terms are ignored. Since the e'^ is defined along the direction of R4, the z coordinate 
corresponds to the shared Hydrogen stretching vibration. 

For the other 6 DOFs of the (0H)2 fragment, the 6D KEO is written in terms of 
the Gvib matrix elements Eq. (1751) . Here we take the angular momentum of the (0H)2 
fragment as total angular momentum since the angular momentum of the shared Hydrogen 
is negligible. Combining both parts we finally obtain the 9D KEO where we introduce 
the new coordinates Ui = cos 6i {i = 1, 2). For the simulation an additional normalisation 
transform was performed to reduce the numerical effort. This gives the following 9D KEO: 

T = Ti + T2 + T3 (30) 



with 



^ 2^1 dRl 2^2 dRl 2fi4dRl ^ ' 
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^ i 2fi,Rf 1 - ^ 2fi4Rl 1 - / dip 



«=1,2 



1 ui U2 d d 

H cos 09 — 

l^iRl a/1 - ul a/1 - ul dip dip 



1 U2 f d . ~ d d r ~ . d \ 



n — ^ o I ^siny^-y/l -Mf— + —Jl-ufsinip— 

^HiRi - ul \oip V dui dui V dip J 

-wj!7r=lU""^v'""=a^^ssvi ■ (32) 

The non-Euclidean normalisation according to the volume element is dj. = 
dRidR2dRidxdydzduidu2dip. The reduced masses for different isotopomers are defined as 
follows: HOHOH^ - /^i = Ai2 = ""^h^o/I^h + ^o), /^s = 2mH(mH + mo) / (3mH + 2mo), 
/i4 = (mn + mo)/2; HODOH" - /^s = 2mD(mH + mo)/(mD + 2mn + 2mo); HOHOD" - 



yUi = mj^mo / {m-ii + mo) and /is and /i4 change correspondingly. For DODOD" the cor- 
responding masses of HOHOH" are modified by replacing mn by m-Q. In the same way 
we get the masses for DOHOD^ by exchanging mn and rriD in HODOH", and similarly 
one can obtain DODOR- from DOHOR-. 



B. Reduced Dimensional KEO for N2H7 

In this section we give an example for a KEO for a system with constraints stressing 
again that our method is invariant for reduced dimensional descriptions. Specifically 
we focus on the protonated ammonia dimer, N2H7', whose IR spectrum in the range < 



1500 cm ^ had been discussed in Ref. 



26| on the basis of a 5D model including only 



the asymmetric stretching vibration of the shared proton. Subsequently, the model was 
extended to account for the two degenerate bendings as well in Ref. |27|. Here we will 
derive the 5D model and comment on its extension. A reasonable reduced description for 
the spectrum which is infiuenced by the shared proton motion should take into account 
the following two conditions: (i) experimental data suggest that the relevant energy range 
is presumably below 1000 cm~^ , but might extend into the < 1500 cm~^ range due to 
combination bands and (ii) symmetry selection rules dominate the anharmonic couplings 
especially in this low-energy range. Therefore, we first assume that the C3 symmetry of 
the N2H6 fragment, i.e., excluding the central proton, will not be broken. Second, the 
length of the N-H covalent bonds shall be fixed. Notice that this neglects the bending 
motion of the shared proton perpendicular to the line connecting the two nitrogen atoms. 
This is justified by the different symmetry and the lower anharmonicity of these two 
modes as compared to the proton motion along the N-N line. These constraints leave five 
internal coordinates to describe the system as shown in Fig. [21 i.e., the shared proton 
stretching with respect to the center of mass of the rest N2H6 fragment, the relative 
motion of the centres of mass of the two ammonia, i?, the umbrella type motion of the 
two ammonia, 9i and ^2, and the rotation (torsion) of the NH3 fragments with respect to 
each other, ip. 

For these five active coordinates the classical kinetic energy with all other internal 



coordinates frozen can be obtained by exploiting Eq. (fT5|) repeatedly 

rp _ ^(NaHe) ^.(NaHe) 
^ ~ ^(H) ^ ^NaHe 

_ ^(NaHs) , ^(NHa) , ^(NH?,) ^(NH?,) 

~ ^(H) "^^NHs "^^NHJ, "^^(NHg) 

_ ^(NaHs) ^(NHJ,) /^(Hg) ^(H3)\ /^(H?,) ^(Hi)\ .ooX 

~ ^(H) ^(NH3) l^^Hs "^^(N) y l^^HJ, "^^(N') l'^'^^' 

where (AB) is the center of mass of AB. It is straightforward to obtain each term in 
above equation 



(N2H6) _ 1 ,, ^2 
2 ' 
1 



-^(H) ~ 2^^^ 



^(NHs) _ -L p2 
^ (NH3) - 

(H3) 3 /d(i?osin6'j)V 3 • /)\2-2 

^H3 = 2"^H ( Ji J + 2"^H(i?osm^i) V'i 



.(H3) _ 1 ,. diRoCosOi 



2 



^« -^<-H,^^i^j . (34) 

where /Xp = 2mH(3mH + mN)/(7mH + 2mN), /i/j = ^(3mH + mt^) and /i(N-3H) = 
Sm-nrri'^ / [SniYi + mN). The orientation angle of each individual ammonia is denoted by 
ipi^2 and only the difference between them is the torsion shown in Fig. [2] and i?o is the free 
N-H covalent bond length. Thus the final kinetic energy is 



T = ]-^,z' + ^^aR' + j^"^^^^ Rl ((^^,ln'e, + ^.^sin^^^ 
2 2(3mH + mNj V 

+ ^-mnRlY,{Oho^'()^ + ^h^^%) 



i=l,2 

1 .0 1 1 . .AO 1 



2 



+ ^J,ot(^l)0? + ^/rot(^2)0^, (35) 

where /vib(^) = /o(cos^6' + mNsin^6')/(3mH + m^), Imti^) = losin^O and Iq = 3mH-Ro- 
The last two terms in above equation can be rewritten and after separation of the global 
rotation we finally simplify the kinetic energy as 

T = ^fipz^ + ^finR^ + l^Mn^ + ^Aor(^i, ^2)0', (36) 

i=l,2 

where the reduced moment of inertia for the torsion is /tor (^1,^2) = 

/rot(^l)/rot(^2)/(/rot(^^l) + /rot(^^2)). 

So far we did not consider the Coriolis type couplings with total angular velocity. 
However, due to the C3 symmetry, the total angular momentum equals to zero is equivalent 



to total angular velocity equals to zero for this specific case since none of the adopted 
coordinates contributes to total angular momentum. Therefore the Coriolis type couplings 
are zero. 

According to the general procedure detailed in Section Hi Al we can obtain the following 
quantum KEO {ip = (p2 — ^i) 

1 - 1 - 1 - 1 - P^P-^ 

+ lPiMi0i)Po^ + lpSMi02)Pe,. (37) 

Inspecting the adopted coordinates we notice that one of them is Cartesian and the other 
four are spherical coordinates defined in MRF. Therefore, according to Section III Dl the 
HCMOs are given by 

Pi = Pz 

Pr - Pb-^ 

Pl^ = Pe^-incoiOj,] = 1,2 

Pi = P^- (38) 
Hence the final KEO and the Euclidean normalisation condition are given as follows 

d"^ _h'^ ^ 1 d sine, d 
~2IU0MW '^h'a 
dr = R"^ sin 9i sin 92dzdRd9id92dip. (39) 

In Ref. 



27| we have extended this model by including the shared proton bending 
coordinates x and y as shown in Fig. [2l In terms of the KEO this is straightforward since 
in Eq. flMj) we simply have to set 

Tlu^'' = ll^pi^' + + f) (40) 

which is carried through the derivation to give a corresponding contribution to Eq. fl39l) . 
Compared to the experiment this gives a rather good description of the bending funda- 
mental transitions as well as an improvement concerning the shared proton fundamental 
stretching transition. 



IV. SUMMARY 



In summary we have presented a new method for obtaining kinetic energy operators 
which involves a straightforward quantisation of the classical Lagrange/Hamilton function. 
It is based on the notion of hermitian conjugate momentum operators whose derivation 
presents the major effort for practical applications. An important point of our approach is 
its validity for systems with constraints. This makes it particularly attractive for studying 
larger systems in reduced dimensionality. 

We have applied our approach to the case of a general system described by polyspherical 
coordinates and showed the equivalence of the kinetic energy operator with that obtained 
by the method of Gatti and coworkers. Afterwards we discussed two specific applications 
to charged cluster systems having strong hydrogen bonds. The strong coupling between 
different coordinates as well their structural floppiness makes it necessary to perform high 



dimensional quantum dynamics simulations as done in Refs. 
the kinetic energy operators derived in Section lllli 



25 



26 



27| on the basis of 



V. APPENDIX 



A. Hermitian Conjugate Momentum Operators for Spherical Coordinates 

The derivation of the expressions for the HCMOs in LRF is the most tedious part of 
the present method and it becomes even more involved if coordinates defined in MRFs 
are used. In this Appendix we will first summarise the relevant relations between the 
LRF and MRFs and subsequently turn to the derivation of HCMOs. 

Consider the LRF defined by three orthogonal unit vectors {e^,, 6^^,62} and a MRF 
defined by three orthogonal unit vectors {e^;', e^/, e^/}. The orientation angles of e^/ in 
the LRF are ('(9,0). To obtain the connection between the LRF and the MRF we first 
apply two excessive rotations Uy{'&) and Uz{4') to the LRF, where Vy{0) means rotating 
around By and Uz{4>) means rotating around e^. The matrix representation of them 



are 



^ cos?? sin-i? ^ 



U,(^) 



1 
— sin cos 



cos (j) — sin (f) 

U^(0) = sin0 COS0 . (41) 

V » "V 

The new reference frame generated by applying Uy('(9) and '[Jzi'P) to the LRF has the 
same z axis as the MRF, i.e., the only difference between the two reference frames is 
just a rotation by x around e^/. Therefore, the MRF can be obtained by applying three 
successive rotations \Jy{;d), \Jz{.'P) and '^z'{x) to the LRF, namely 



U.'(x)U,(0)U,(79)e, 



a = x,y,z, 



(42) 



where \Jz'{x) means rotating x around e^'. The matrix representation for \Jz>{x) in the 
MRF is the same with Uzix) m the LRF. Since the third rotation \Jz'{x) in Eq- (1421) 
has no effects on e^/ we actually have e^/ = U z{(t>)Uyi^^)e;^. Exploiting the rules between 
vector and operator transformations we can obtain the expression for \Jz' (x) in the LRF 



Vz'ix) = [U,(0)U,(^)] Vzix) [Vz{<jy)Vy{^)]-' . 



Consequently Eq. fl42l) can be rewritten as 



U^(0)Uy(^?)U^(x)e„ a = x,y,z. 



(43) 



(44) 



Eq. (jH]) tells us that an equivalent way to obtain the MRF is to apply three successive 
rotations U^(x), Uy(?9) and lJz{4>) to the LRF. 

Now let us consider vectors Rj characterised by spherical coordinates {Rj,9j,ipj) in 
the MRF (j = 1 — A^). With Eq. (jUj) and the rule between basis vectors and components 
transformations we can obtain the Cartesian components of these vectors in LRF 



R 



'jy 



U,(0)U,(^9)U,(x) 



^ Rjx' ^ 

Rjyl 



(45) 



where Rj^ and Rj^i' are Cartesian components of Rj in the LRF and the MRF, respec- 
tively. We can rewrite Eq. (1451) in more formally as 



R, = R,Vz{4>)Vy{^)Vz{x)Vz{v,)U: 



(46) 



Eq. (H6|) is a vector equation therefore it is also valid for an arbitrary reference frame. 

Now we turn to the main task, that is, deriving the Hermitian conjugate of momentum 
operators. This can be done step by step according to Eq. (jS]) with the help of the 
coordinate transformation Eq. (jlSD- From Eq. we can see that {6j, ipj) only appear in 
the Cartesian components of Rj. This greatly simplifies the expression for the Hermitian 
conjugates of momentum operators as follows 



a=x,y,z 



a='x,y,z 



d dRja 
dRja dRj 

d dRja 
dRja dOj 



a=x,y,z ^ ■' I J / 

Here the three Euler angles in Eq. (H5|) are purely parameters. To calculate the partial 
derivatives we first give some useful relations following from the orthogonality of the 
transformation Eq. (145 p 

dRja p p _ ^ dRjp, 

- 2^ m;;,^^^' ' - -dR~a 



2 



dRja \ \^ { dRja 



dRja dRjfji 



(4J 



dRjfji dRja 

where a = x,y,z and /5' = x',y',z'. Note here all the derivatives {dRj^i / dRja} are 
parameters only depending on the three Euler angles. 



Let us first derive the expression for Pt . It is straightforward to get the following 



derivatives of MRF Cartesian coordinates 

dR-jx' ■ ri Rjx' 

= sm &j COS (pj = 
dRj Rj 

dRjyl . . Rjyl 

= sm tij sm (fj = 
dRj Rj 

where i?| = + + = + + i?!^,. Based on Eq. (jH]) we can obtain the 



derivatives of LRF Cartesian coordinates with respect to Rj and derive the final result 

E/ d dRja \° _ sr^ I 9 Rji3' dRjc 
\dR~~dR~) ~ ^ I aR~^~R~aR^/ 




1 OR,. HR% + R% + R%r\ dR 



Rj dRja dRja j ORjfS' 

1 dRjfji RjfjiRja \ dRja 



Rj dRja Rj J dRjfji 

Rja 

J-Lj JXj JXj 

where Eq. ( 148 p has been used in the fourth step. Finally we can get the HCMO associated 
with Rj as follows 

Next we will derive the expression for following the same procedure. The derivatives 
of MRF Cartesian coordinates read 

dRjx' „ „ 

^— — = Rj cos 6j cos (pj = Rjx' cot 6j 
dOj 

dRjy' 

— --j — = Rj cos 6j sin ipj = Rjyi cot Oj 
dOj 

—77— = -Rjsm6j = (52) 

dOj ^ ^ coiOj ^ ' 

Based on Eq. ( 152|) we can obtain the derivatives of Cartesian coordinates in LRF with 
respect to Qj 

Furthermore, we obtain 

d dRja \ „ dRja dRjx' ^ dRja dRjy' 1 dRja dRj^' 



dRja dOj J dRjx' dRja ^Rjy' dRja cot 9j dRjz' dRjc 



^ dRja ^ ^ dRja ^ Rjz' dRja \ dcotOj 
dRjx' "'^ dRjy' cot^ 9j dRjz' ) dRja 

cot Q ^ dRja dRjz' 

sin 9j cos 9j dRjz' dRja 

( / 9 ^ X dRja \ d cot Qj , , 

+ \R,a + (tan^ - 1)^..'^ j -Q^- (54) 



Since cot6'j = Rjz' / J R^j^, + R^^y, we can get the following derivatives 



dcotOj Rjx'Rjz' 

dRjx' fW2 I d2 

d cot 6j RjyiRjz' 

dRnii' 



d cot 9j 1 



(55) 



Based on Eq. ( |55ll we have 



^ ^ R,.R,. + R]^]. (56) 



'p2 _|_ p2 \ " dRja 



Combining Eq. fl54p and Eq. fl56p we can now derive 



E 



9 (9i? 



'J a 



dRja dRjz' 



3 cot 6'j 



sin dj cos 6'j 



+ — 3 (^—Rjz'Rj + R?jRjz' + (tan^ 0j — \){—R?j^iRjz' + R^Rjz')) 

^JR'jx' + ^jy' 

= 3 cot + (tan^ - 1) cot 9^ 

sm 9j cos 

= cot6'j (57) 

The Hermitian conjugate of each Pg. associated with the corresponding MRF polar angle 
6j can be finally expressed as 

moots, = --i^|-=in9, (58) 
Finally we come to P^^. The procedure is the same, however, it is much simpler as 



compared to Pj. . Again we start with the derivatives of MRF Cartesian coordinates 

dRjx' _ 



dRjz' 

to obtain the final 



—Rj sin 9j sin Lpj = —Rjy' 
Rj sin 6j cos (fj = Rj^i 

(59) 




(60) 



dRja dRjyi dRja dRjx' 
The momentum operator associated with each orientation angle (pj is Hermitian according 



to Eq. (1601) . namely 



Let us recall the above detailed procedure. The only condition we need is that the 
transformation matrix between Cartesian coordinates in the LRF and those in a MRF 
is orthogonal. First, if we set all the three Euler angles equal to zero we can obtain the 
HCMOs associated with spherical coordinates in the LRF. The final results are the same 
with Eq. flHT]) . Eq. (1581) . and Eq. (1611) since a unit matrix is also an orthogonal matrix. Sec- 
ond, we can use more rotations to define more MRFs and the HCOMs associated with the 
spherical coordinates in each different MRF obey Eq. (1511) . Eq. ( l58l) . and Eq. (1611) . As an 
important consequence the Hermitian conjugates of all the momentum operators associ- 
ated with real bond lengths, bond angles and dihedral angles obey Eq. (ISTI) . Eq. (158!) . and 
Eq. (|6T]) . respectively, irrespective of how complicated the molecule might be. Similarly, 
the result holds for any reference frame provided there exists an orthogonal transforma- 
tion to transform it to the LRF. Based on the above conclusions we can see the HCMOs 
associated with the three Euler angles obey the same relations 

ih d 





= p^ 


— ih cot -i? 




= h 


= -^^^ 


X 


= Px 


= -in— 
ox 



sin -(9 d-& 



sin ^ 



(62) 

Finally, we can obtain the HCMOs Pj^, associated with MRF Cartesian coordinates 
following the same procedure. Since dRja/ dRja' is just a parameter which does not 



depends on Rj^, Rjy or Rj^ we can easily get 

(63) 

where a' = x',y',z'. Apparently, the momenta associated with MRF Cartesian coordi- 
nates are Hermitian. 



B. KEO in Terms of Polyspherical Coordinates 

In order to connect our approach to existing one we consider the KEO in terms of 
so called polyspherical coordinates defined in the LRF and a MRF. This will lead us to 
expression which have been reported in Ref. |l7|. Consider a molecular system composed 
of + 1 atoms. After separating the total centre of mass motion we can describe it with 
A^ vectors -Ri, -R2, ■ ■ ■ , Rn- The three Euler angles {{}, 0, x) are chosen in such a way that 
the spherical coordinates of in the LRF are {Rn, ■}},(/)) and the spherical coordinates 
of Ri in the MRF are {Ri, 61, (fi = 0). That is to say the e^,/ axis of the MRF is defined 
to be along the direction of Rn- The rest N — 2 vectors are characterised by spherical 
coordinates {Rj, 9j, ipj} in the MRF (j = 2, ■ ■ ■ , A^ — 1). 

Using Eq. ( l20i) with j = 1, - ■ ■ ,N and 9n = fN = V^i = 0, we can derive the LRF 
components of the A^ velocity vectors 

R, = i?,i9U,(</.)U;(79)U,(x)U,((^,)U,(^,-)e. 

+ RjVM^yi^)^z{x)^z{v,)Vy{e,)e, 
+ i?,-^,U,(0)U,(79)U,(x)U,(<^,)U;(^,)e, 

+ P,-v9,U,(0)U,(^)U,(x)U;(v9,)U,(0,)e., (64) 

where U' is the first order derivative of the corresponding rotational transformation ma- 
trix. The above equations are still vector equations, however, they are only valid in the 
LRF since we have used the relation = 0. To make them valid in an arbitrary reference 
frame we only need to add the corresponding terms containing to the right hand side. 

Before proceeding we will discuss the orthogonality of the terms appearing in Eq. (1641) 
. The last three terms in Eq. fl64l) are just the spherical velocity components measured 
in the MRF while the first three terms are components of the velocity caused by the 



non-inertial MRF. Each term in Eq. (IMI) can be expressed in a formally simple way in 
terms of angular velocity vectors 

itj = de^ X Rj + X Rj + x Rj 

+ RjRj/Rj + OjBQ, X i?j + (fjB^. X (65) 

where is the direction of the angular velocity {}, and similar for e^, e^, e^,, and Ci^^.. 
With the help of Eq. (1651) one can immediately see that the fourth term is parallel to 
Rj while the other terms are perpendicular to Rj. Recalling the velocity in spherical 
coordinates we know that the last three terms are orthogonal to each other. Therefore, in 
Eq. (l64l) (or Eq. ( l65i) ), the fourth term (namely the Rj term) is orthogonal to the other 
terms and the last three terms are orthogonal to each other. 

Having at hand all the velocities it is quite straightforward to write the classical kinetic 



energy 



28l | according to 

1 ^ 

T = -J2 (66) 



where the matrix {fijk} combines the reduced masses and the transformation matrix 



281]. The final results in LRF can be written in a 



between adopted and Jacobi vectors 
symmetric form as in Eq. ([T]) 

T = WmQ, Q = I I , (67) 

^ \ Qvib 

where the coordinates are separated as rotational, Eq. (|22l) . and vibrational DOFs 

Qlib = [Ri ■■■ Rn 0, ■■■ d^_, ^2 ■■■ V>N-i ) • (68) 

The next step is to calculate the inverse matrix of M to generate the quantum KEO. 
Aiming at a separation of the rotational and vibrational DOFs, we divide the matrix M 
into the following four blocks 

/ Mrot Mcor \ 

M = , (69) 

where the subscripts "rot" , "cor" , and "vib" correspond to rotational, Coriolis, and vibra- 
tional terms, respectively. Therefore, Mrot, Mcor, and Mvib are 3 x 3, 3 x (3A^ — 3) and 
(3A^ — 3) X (3A^ — 3) dimensional matrices, respectively. Suppose the inverse matrix of M 
is divided in the same spirit as 

M-^ = I I , (70) 



we can express the quantum KEO as a sum of rotational, vibrational, and Coriolis terms 



T — Trot + T^ih + Tcor 

2 



+ I1.C.). (71) 



According to Section III El the generalised momentum vector Prot associated with the 
rotational DOFs is just the total angular momentum vector J. If we are interested in the 
rotational DOFs only, we obtain upon setting the vibrational momenta equal to zero: 

Trot = 2-^rot^rot-Prot = '^J Grot«^ ; (72) 

where the components of J should be along the directions of e^, and e^^ according to 
Eq. (|27|) . In passing we note that one can also transform the components to Cartesian 
ones according to Eq. (!28|) . On the other hand, if we are interested in the vibrational 
DOFs, the KEO for the total angular momentum J = reads 

Tvih = -Pyih^vihPvih- (73) 

Finally, we give the details on how to calculate the different blocks of the inverse 
matrix of M. One can see that the matrix M can be congruently block-diagonalised in 
the following way 

' 1 -McorM;.^ 
1 




Cor 





Calculating the inverse of both sides of Eq. ( 17il) leads to 

-1 



Mcor 


I ( 








I- 







(74) 



M 



-1 




Mrot Mcor 

Mrot - McorM^i^Mf,^^ 
-M-^mL 1 i I 



-1 




(75) 



^ib^^-^Cor / \ " ^^-"vib 

Calculating the matrix product in above equation leads to the final results 

Grot = (Mrot - McorM^i^M;.^^' 
Gcor = — GrotMcorM^iJ, 

Gvib = M;^^ + M;i^M{,„,Gcor. (76) 



The matrix G is essentially the same as the one reported in Ref. [l7j. There, however, 
the essential point has been to express the MRF Cartesian components of the angular 
momentum associated with each vector {Ri} in terms of J and {Pe-,P<^.}. 

For the special case of Jacobi vectors, M^ib is diagonal which greatly simplifies the 
calculation of Eq. (176 p . It is straightforward to obtain the diagonal elements according to 
Eq. dMl) 



(Mvib)^^,,^. =/i,i?|W%. (77) 

The inverse matrix of Mvib can be obtained quite easily. We only need to calculate the 
inverse of a matrix of dimension 3x3 as well as some matrix products to get the final 
results. Here we give the matrix elements for Gvib (Hermitian) which will be used in 
Section Hm 

Gr^Rj = — , Gr-q. = 0, Gr-^^ = 0, 

/ij 

6ij cos((^i - ipj) 



G, 



cot 9j sin(y9j — y^j) — cot 9i sin ipi 



_ 



fijRj sin^ 9j fiiRf sin^ 9i 
^ cot 9i cot 9j cos{ipi — (Pj) + cot^ 6*1 — cot 9i (cos ipi cot 9i + cos (fj cot 9j) 
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Figure Captions 



FIG. 1: The four Jacobi vectors of the H3O2 anion which are used to define the 9D KEO. 



FIG. 2: Definition of active coordinates of the reduced N2H7 model. 




Figure 1, Yang and Kiihn 




Figure 2, Yang and Kiihn 



